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ABSTRACT 

The discovery by Swift that a good fraction of Gamma Ray Bursts (GRBs) have a 
slowly decaying X-ray afterglow phase led to the suggestion that energy injection into 
the blast wave takes place several hundred seconds after the burst. This implies that 
right after the burst the kinetic energy of the blast wave was very low and in turn 
the efficiency of production of 7-rays during the burst was extremely high, rendering 
the internal shocks model unlikely. We re-examine the estimates of kinetic energy 
in GRB afterglows and show that the efficiency of converting the kinetic energy into 
7— rays is moderate and does not challenge the standard internal shock model. We also 
examine several models, including in particular energy injection, suggested to interpret 
this slow decay phase. We show that with proper parameters, all these models give 
rise to a slow decline lasting several hours. However, even those models that fit all 
X-ray observations, and in particular the energy injection model, cannot account self- 
consistcntly for both the X-ray and the optical afterglows of well monitored GRBs 
such as GRB 050319 and GRB 050401. We speculate about a possible alternative 
resolution of this puzzle. 

Key words: Gamma Rays: bursts— ISM: jets and outflows-radiation mechanisms: 
nonthermal— X-rays: general 



1 INTRODUCTION 



The X-ray telescope (XRT) on board Swift has provided high 
quality early X-ray afterglow light curves of many Gamma- 
ray Bursts (GRBs). One of the most remarkable and un- 
expected features discovered by Swift was that many of 
these X-ray afterglow light curves are distinguished by a 
slow decline — The flux F decreases with observer's time t 
as F oc ^I * -0 - 8 ^ lasting from a few hundred seconds to few 
hours (Nousek et al. 2005, Campana et al. 2005; Vaughan 
et al. 2005; Cusumano et al. 2005; de Pasquale et al. 2005). 
Such a phase is unexpected in the standard fireball model. 
A simple explanation is that the slow decline arises due to 
a significant energy injection (Zhang et al. 2005; Nousek et 
al. 2005; Panaitescu et al. 2006, Granot & Kumar 2006), 
as suggested previously (For baryon-rich injection: Rees & 
Meszaros 1998; Panaitescu et al. 1998; Kumar & Piran 2000; 
Sari & Meszaros 2000; Zhang & Meszaros 2002; Granot, 
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Nakar & Piran 2003. For Poynting flux dominated injection 1 : 
Dai & Lu 1998a; Zhang & Meszaros 2001; Dai 2004). It 
has been argued that consequently the resulted GRB effi- 
ciency, i.e., the ratio of the energy emitted in 7— ray energy 
to the total energy (the sum of the 7— ray energy and the ki- 
netic energy of the ejecta powering the afterglow), should be 
90% or higher. Some extreme assumptions are needed (Be- 
loborodov 2000; Kobayashi & Sari 2001) to reach such a high 
efficiency within the framework of the standard internal- 
shocks model (Paczynski & Xu 1994; Rees & Meszaros 1994; 
Sari & Piran 1997a , 1997b; Kobayashi, Piran & Sari 1997; 
Daigne & Mochkovitch 1998; Piran 1999). 

We re-examine this issue focusing on two critical aspects 
of the analysis. The estimate of the kinetic energy of the 



If the outflow ejected from the central engine after the gamma 
ray burst phase is highly magnetized, at a radius ~ 10 15 cm, 
the MHD condition breaks down. Significant magnetic field dis- 
sipation processes are expected to happen which converts energy 
into radiation. As long as the highly magnetized outflow is steady 
enough, strong and slowly decaying X-ray emission is possible (see 
Fan, Zhang & Proga [2005a] and the references therein). 
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ejecta from the afterglow observations and in particular from 
the X-ray flux and the need of energy injection. We show in 
§2 that even for these Swift GRBs with long duration X- 
ray flattening the 7-ray conversion efficiency is high but not 
unreasonable. 

We then turn to the puzzling slow decline seen in the 
first few hours of the X-ray afterglow. We explore in §3 sev- 
eral models that may give rise to slowly decaying X-ray af- 
terglows: (i) Energy injection, (ii) A small £ e , in which only 
a small fraction, £ e <C f of the electrons are accelerated to 
high energies and contribute to the radiation process, (iii) 
Evolving shock parameters, where the microscopic shock pa- 
rameters e e and/or es (the fraction of shock energy given 
to the magnetic filed) vary in time and are inversely propor- 
tional to the Lorentz factor of the ejecta. (iv) A very low 
variable external density model, in which the number den- 
sity of the medium is not only very low but it also a function 
of the radius, (v) Highly magnetized outflow where flatten- 
ing might arise because of a slow conversion of the magnetic 
energy to kinetic energy of the external matter. We present 
in §3 analytical derivation as well as numerical calculations 
of the expected light curves in all these models except the 
last one. In §4 we compare the models to the observations 
of GRB 050319 and GRB 050401. We summarize our re- 
sults and discuss their implications in §5. We conclude with 
a speculation on the nature of the solution to this puzzle. 



2 IS THERE A GRB EFFICIENCY CRISIS? 

One of the critical factors that characterize the emitting of a 
GRB is the energy conversion efficiency. The 7-ray efficiency 
is defined as: 

e 7 = E 1 /{E 1 + E k ) , (I) 

where E 1 is the isotropic equivalent energy of the 7— ray 
emission and Ek is the isotropic equivalent energy of the 
outflow powering the afterglow. Following the Swift obser- 
vations of flattening in the X-ray afterglow light curve of 
many GRBs, it has been argued that typical values of e 7 
could be as high as 90% or even higher (Zhang et al. 2005; 
Nousek ct al. 2005; Ioka et al. 2005; for the discussion of 
pre-Swift GRBs, see Llod-Ronning & Zhang 2004, hereafter 
LZ04) . This very high efficiency would challenge most 7-rays 
emission models and in particular it challenges the standard 
fireball model that is based on internal shocks. 

These claims arise from revised estimates of the kinetic 
energy immediately following the GRB. Therefore, in or- 
der to explore this issue we re-examine the estimates of the 
kinetic energy from the X-ray observations. As we show be- 
low, at a late afterglow epoch, the X-ray band is above the 
cooling frequency. In this case the X-ray flux is independent 
of the poorly constrained n and the X-ray luminosity is a 
good probe of Ek (Freedman & Waxman 2001, Kumar 2000, 
LZ04). 

In the standard GRB afterglow model (e.g., Sari, Piran 
& Narayan 1998; Piran 1999), the X-ray afterglow is pro- 
duced by a shock propagating into the circum-burst matter. 



The equations that govern the emission of this shock are 
(Yost et al. 2003) 2 

F„, max = 6.6 mJy (^^)Dl 2 28 34 e 1 J 2 _ 2 Ek } 5 3 n /2 , (2) 

, c „ 1n llTi Eil/2 1/2 2 W2/ 1 + 2x1/2.-3/2 /„n 

v m = 7.6 x 10 Hz S k / 53 e B / i _ 2eei _ 1 C p (— — ) t d , (3) 

i ^„inl5tT c-1/2 -3/2 -1 / 1 + 2 n -1/2 .-1/2 1 , .s 

v c = 1.4x10 Hz E k ^ e B /_ 2 n (—3—) t d + y) 2 '( 4 ) 

where z is the redshift, Dl is the corresponding luminos- 
ity distance, p is the power-law index of the shocked elec- 
trons, we use p = 2.3 throughout this work, C p = 13(p — 
2)/[3(p — 1)] and td is the observer's time in unit of days. 
Y — (—1 + \/l + 4w JCJV e e /es)/2 is the Compoton param- 
eter, where r\ = min{l, {f m /i'c) < ' p ~ 2 ^ 2 } (e.g. Sari, Narayan 
& Piran 1996; Wei & Lu 1998, 2000), < r] KN < 1 is a co- 
efficient accounting for the Klcin-Nishina effect, which is -y E 
(the random Lorentz factor of the electron) dependent (see 
Appendix A for detail). Here and throughout this text, the 
convention Q x — Q/1GT has been adopted in cgs units. 

For the typical parameters taken here, u m crosses the 
observer frequency vx ~ 10 17 Hz at td ~ 4x 10~ 4 . It is quite 
reasonable to assume vx > max{f c , z/ m }, and the predicted 
X— ray flux is 

F - F v 1/2 u^- 1) ' 2 v- p/2 

= 3.8 x 10- 4 mJy (^ 2+ * /4 D^ 8 . 3i e% ~_T^\ 

El% 2)/ \l + Y)-¥t 3p)/4 «- x % 2 . (5) 
The flux recorded by XRT is 

T = / F Vx dvx 

Jv X l 

= 1.2 x 10- 12 ergs s" 1 cm" 2 (l^l^^D^ 
4 P - 2 2 )/4 C-i473 2)/4 (l + Y)-Ht Zp)/ \ (6) 

2 To derive these equations, the deceleration of the fireball is gov- 
erned by the energy conservation F 2 Mc 2 = Ek, where M is the 
rest mass of the shocked medium (e.g., Blandford & McKcc 1976; 
Sari et al. 1998; Piran 1999). The distribution of the fresh elec- 
trons accelerated by the shock is assumed to be dn/d-je oc ■y^ p 
for 7 e > 7e,m, where 7 e , m = (m p /m e )[(p - 2)/(p - l)]e e (r - 1), 
governed by the strict shock jump conditions (Blandford & Mc- 
Kee 1976). The other crucial parameter is the cooling Lorentz 
factor 7e jC = 6(1 + z)nm e c/[axrB 2 t(l + Y)], above which the 
energy loss due to the synchrotron/invcrse-Compoton radiation 
is important (Sari et al. 1998; Piran 1999), where ctt is the 
Thompson cross section and B is the magnetic field of the 
shocked medium. v m and v c are the corresponding synchrotron 
radiation frequency of electrons with Lorentz factor 7 e ,m and 
7e,c: respectively. The maximum specific flux is estimated as 
F„, max Rb (1 + z)Mre 3 B/(4:Tvm p m e c 2 Dl ) (Sari ct al. 1998; Wi- 
jers & Galama 1999), where e is the charge of electron. The v c 
and -Fj/.max taken here are comparable with that of most previous 
works (e.g., Granot et al 1999; Wijers & Galama 1999; Panaitescu 
& Kumar 2002; LZ04). The v m is close to that taken in Sari et 
al. (1998), Granot ct al. (1999) and Wijers & Galama (1999), but 
is about 30 — 40 times smaller than that taken in Panaitescu & 
Kumar (2002) and LZ04. Such a large divergency may arise if one 
ignores the term (p — 2)/(p — 1) when evaluating 7 e ,m- 
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where vxi = 0.2 keV and vxi = 10 keV. This equation is 
now inverted to obtain Ek from the observed flux. 

In some special cases, v m < vx < v c , the flux recorded 
by XRT should be 

T = 1.5 x 10- 11 ergs s" 1 cm" 2 {±±1)^/^1% , 34 

,(P+1)/4,P-1 1/2 p (p+3)/4,3(l-p)/4 (7 , 
fc B,-2 fc e,-l' l ^£,53 L d ' \'J 

2.1 The efficiency of the pre-Swtft GRBs 

With equation (6), the corresponding X-ray luminosity at 
td — 0.4 (~ lOh, to compare with the results of LZ04) is 

L x = AtxD\T/{1 + z) 

i i „ i n46 -1 -2 / 1 + Zn(p-2)/4 

= 1.1 x 10 ergs s cm ( — - — y p " 

4 p : 2 2 )/4 <-- 1 1 (i+^)- 1 4^ 2)/4 1 (8) 

which in turn yields 

E k = 9.2 xlO 52 ergs i?L^ ( 4 p + 2) (i±i)( 2 - p )/( p + 2 ) 

£ -(p-2)/(p+2) £ 4(l-p)/(p+2)^ 1 + Y^i/(P+2) (g) 

where 7? ~ [t(10h)/TQ ] 17e '^ 16 is a factor accounting for the 
energy radiative loss during the first 10 hours following the 
prompt gamma-ray emission phase (Sari 1997; LZ04), T90 
is the duration of the GRB. The numerical factor of our 
equation (9) is larger than that of equation (7) of LZ04 by 
a factor of 9.2(1 + Y) i/(p+2) due to the facts that (1) The 
v m taken here, which matches the numerical result better 
(one can verify this with a simple code to calculate the dy- 
namical evolution as well as v m numerically), is about one 
and half orders smaller than that taken in LZ04. (2) The 
inverse Compton effect has been taken into account. Similar 
conclusions have been reached by Granot, Konigl & Piran 
(2006). However, it is not easy to estimate Y since it de- 
pends on es sensitively (see Appendix B for discussion). 
One good way to estimate the GRB efficiency may be to 
take Y ~ and (1 + Y) ~ (e e /e B ) 1/2 , respectively. In both 
cases, our estimates of e 7 (Table 1) are significantly lower 
than those of LZ04. Smaller e 7 may be possible in view 
of that both e e and es might be significantly lower than 
the standard parameters taken here (Panaitescu & Kumar 
2002). We suggest that the typical GRB efficiency of these 
pre-Swift bursts is ~ 0.1 (see Table 1 for detail). Such values 
are well understood within the internal shock model. 

Additional support for this conclusion arises from late 
energy estimates. Berger, Kulkarni & Frail (2004) used the 
late time radio observation to estimate the kinetic energy 
at this stage. The find high energies and correspondingly 
low 7-ray efficiency. For GRB 970508 and GRB 970803, the 

3 While our results are very close to the recent calculations of 
Granot, Konigl & Piran (2006), they also show that the estimates 
of -Bfc arc very sensitive to the exact expressions used for u c , 
u m , and -F^max- Similar conclusion can be drawn by comparing 
previous results of Granot et al (1999), Wijers & Galama (1999), 
Freedman & Waxman (2001), Panaitescu & Kumar (2002) and 
LZ04. Therefore, an alternative explanation for the apparent high 
efficiencies is that the blast wave energy estimates using Lx arc 
simply inaccurate. 
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Figure 1. X-ray (0.2-10 keV) afterglow light curves: Analyti- 
cal (dashed line) lightcurve, and numerical (solid line) when In- 
verse Compoton effect has been ignored. The divergence is about 
a factor of 2. Numerical estimates when the inverse Compoton 
effect has been taken into account with (dotted line) and with- 
out (dashed - doted line) a Klein-Nishina correction. Clearly, the 
Klein-Nishina correction is unimportant for the fiducial parame- 
ters listed in the figure. 



efficiencies are 0.03 and 0.2 respectively, which coincide with 
our estimates (see Table 1). 

The coefficient of our equation (9) are very different 
from that of equation (7) of LZ04. Below we check its valid- 
ity numerically. The code used here has already been used 
in Zhang et al. (2005) and has been tested by J. Dyks in- 
dependently (Dyks, Zhang & Fan 2005). Here we just de- 
scribe briefly the technical treatment. The dynamical evo- 
lution of the outflow is calculated with the formulae pre- 
sented in Huang et al. (2000) , which are able to describe the 
dynamical evolution of the outflow in both the relativistic 
and the non-relativistic phases. The electron energy distri- 
bution is calculated by solving the continuity equation with 
the power- law source function Q — K^~ p , normalized by 
a local injection rate. The cooling of the electrons due to 
both synchrotron and inverse Compton (Moderski, Sikora 
& Bulik 2000) has been taken into account. 

Fig. 1 depicts the numerical results. One can see that 
the numerical results match the analytical ones to within 
a factor of 2. We therefore conclude that equation (6) and 
equation (9) are reasonable approximations to the full solu- 
tion of the problem. 



2.2 The GRB efficiency of Swift GRBs with X-ray 
flattening 

Early flattening is evident for a good fraction of the X-ray af- 
terglow light curves recorded by the Swift XRT. Determina- 
tion of the GRB efficiency of these GRBs is quite challenging 
since, as we see in §4 the underlying physical process that 
controls the slow decline is unclear. A common interpreta- 
tion for this flat decay is energy injection, which essentially 
increases the required initial GRB efficiency. In spite of the 
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Table 1. GRB energies and efficiencies, Lx used in equation (9) and -E 7 are all taken from LZ04. The numerical values quoted in 
parentheses are for (1 + Y) ~ (fe/eg) 1 ' 2 . 



GRB 


_E 7 /10°^ergs 


B fc /10 52 ergs 


efficiency e 7 


970228 


1.42 


17.5 (47.5) 


0.08 (0.03) 


970508 


0.55 


9.1 (24.8) 


0.06 (0.02) 


970828 


21.98 


37.4 (101.5) 


0.37 (0.18) 


97f2f4 


21.05 


78.0 (212) 


0.21 (0.09) 


980613 


0.54 


11.2 (30.5) 


0.05 (0.02) 


980703 


6.01 


22.2 (60.2) 


0.21 (0.09) 


990123 


143.8 


186.6 (507) 


0.43 (0.22) 


990510 


17.6 


121.1 (329) 


0.13 (0.05) 


990705 


25.6 


3.1 (8.5) 


0.89 (0.75) 


991216 


53.5 


337.1 (916) 


0.14 (0.06) 


000216 


16.9 


4.6 (12.5) 


0.78 (0.58) 


000926 


27.97 


91.7 (249.3) 


0.23 (0.1) 


010222 


85.78 


209.7 (569.8) 


0.29 (0.13) 


011211 


6.72 


12.1 (33) 


0.36 (0.17) 


020405 


7.2 


42.3 (115) 


0.15 (0.06) 


020813 


77.5 


203.9 (554) 


0.28 (0.12) 


021004 


5.56 


76.8 (208.8) 


0.07 (0.03) 


XRF 020903 


0.0011 


0.09 (0.25) 


0.01 (0.004) 



uncertainties concerning the applicability of this model we 
consider its implication to the efficiency. 

The energy injection is characterized by a factor / such 
that fEk (/ ~ a few— ten, in the following discussion, we 
take / = 5) is the energy injected into the fireball (Zhang et 
al. 2005). The initial GRB efficiency should be 



e 7 ee E 1 /{E 1 + E k ) = /e 7 /[l + (/ - l)e 7 ] (10) 



where e 7 ee E 1 /(E 1 + fEk) is the GRB efficiency derived 
at td ~ 0.4. LZ04 find that e 7 > 0.4, and therefore, e 7 > 
0.8, which is too high within the framework of the standard 
(internal-shocks) fireball model. However, as shown in §2.1, 
e 7 presented in LZ04 has been overestimated significantly. 
We suggest that e 7 ~ 0.1, therefore even when correcting 
for the additional energy e 7 ~ 0.3, which is still consistent 
with this model. 

As an example we consider the 7-ray efficiency of GRB 
050319. Both the optical (Mason et al. 2005) and the X-ray 
(Cusumano et al. 2005) light curves are well recorded for this 
burst and can be used to constrain the efficiency (see §4.1 
for a detailed discussion). (1) The time averaged optical- 
to- X-ray spectrum (i ~ 200 — 900 s) is a single power law 
with an index (5 = —0.8 (Mason et al. 2005). This implies 
that v m (t ~ 100 s) < v R = 4.3 x 10 14 Hz. (2) The very 
early R-band observation suggests that F„ ima x(t ~ 100 s) ~ 
1 mjy (assuming that energy injection takes place at t > 
400 s). (3) v c > v x ~ 10 17 Hz holds up to t ~ 10 6 s, as 
suggested by the XRT spectrum. We have (see equations 
31-33) e e ~4x 10~ 2 , e B ~4x 10~ 5 , and E h ~ 1.3 x 10 54 
ergs (the energy carried by the initial outflow). With K- 
correction, the isotropic energy of the 7— ray emission of 
GRB 050319 is £ 7 ~ 1.2 x 10 53 ergs (Nousek et al. 2005), 
so e 7 = E 7 / (E~/ + Ek) ~ 0.1. It is sufficiently low to be well 
understood within the standard fireball model. 



3 MODELS FOR A SLOWLY DECAYING 
X-RAY AFTERGLOW 

We turn now to explore (both analytically and numerically) 
models that can give rise to a slowly decaying X-ray after- 
glow phase. The models we discuss include: (i) Energy in- 
jection, (ii) A small Q e . (iii) Evolving shock parameters, (iv) 
A very low nonconstant circum-burst density. We also ex- 
amine the possibility of the X-ray flattening is attributed to 
a highly magnetized outflow. In the numerical calculations 
that we present the parameters are chosen to reproduce the 
XRT light curve of GRB 050319 (for t > 380 s). We also 
present the corresponding R-band light curve. 

3.1 Energy injection 

In the standard fireball model, the fireball that is sweeping 
the circum-burst matter decelerates and its bulk Lorentz 
factor evolves with the time as V oc i~ 3,/8 . With continu- 
ous significant energy injection, the fireball decelerates more 
slowly and slowly decaying multi-wavelength afterglows are 
expected. This model has been analytically investigated by 
many authors (Sari & Meszaros 2000, Zhang et al. 2005; 
Nousek et al. 2005; Panaitescu et al. 2006; Granot & Kumar 
2006). As shown in Zhang et al. (2005), for dE inj /dt oc r q 
we find v m oc t^ 2 ^' 2 , v c oc t (q - 2)/2 , and F„, max oc t 1 " 9 . In 
this subsection, we take q = 0.5 and find: 

{ t (8-7 q) /6 „ f 0.7 5j for y<Vc< Vm . 

t (2-3,)/4 _ t l/ 8j for 
V c < V < V m \ 
t (8-6,)/6 „ t 0.92 j for Vc<v< Vm . (11) 

t [(6-2p)-(p+3),]/4 _ t -0.S!l j for Vm<v< „ c . 
t [(4-2p)-(p+ a ),]/4 _ t -0.88 j for y > max{ ^ c) Um} 

Following Zhang et al. (2005), we consider an energy 
injection rate of the form (1 + z)dEi„j jdt = Ac 2 (t/to)^ 9 
for to < t < t e , where A is a constant. With the energy 
injection, the equation (8) of Huang et al. (2000) should be 
replaced by (see also Wei, Yan & Fan 2005) 
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Figure 2. The X-ray (0.2-10 keV) afterglow light curve and the 
R-band light curve for the energy injection model. 



Figure 3. The X-ray (0.2-10 keV) afterglow light curve and the 
R-band light curve for the small f e model. 



dr 



(1 - V 2 )dm + A(t/t Q y q [dt/{\ + z)] 
M e , + em + 2(1 - e)Tm 



(12) 



where M e j is the rest mass of the initial GRB ejecta, m is 
the mass of the medium swept by the GRB ejecta, which 
is governed by dm = A-nR 2 nm p dR, m p is the rest mass of 
proton, dR = F(r + yT 2 - l)cdt/(l + z), e = i]e e is the radi- 
ation efficiency. Our numerical results, the R-band emission 
and the 0.2-10 keV emission, are shown in Fig. 2. 



(see also Zhang et al. 2005), which is inconsistent with most 
XRT observations (Nousek et al. 2005). In addition, in this 
model, the spectral index of the XRT afterglows in the slow 
decline phase is —1/2. It is much harder than that of most 
Swift X-ray afterglows (see Table 1 of Nousek et al. 2005). 
The Swift observations therefore provide us robust evidences 
of that significant part of, rather than a small fraction of 
electrons, have been accelerated in the shock front. 



3.2 Small ( e 

In the standard afterglow model, it is assumed that a frac- 
tion e e of the shock energy is given to all the fresh electrons 
that are swept by the shock front. However, it is possible 
that only a fraction £ e of fresh electrons has been acceler- 
ated, as suggested by Papathanassiou & Meszaros (1996). 
With this correction, equations (2) and (3) take the form 



6.6 mjy ( e 



A + z 



D 



1/2 



Z,,28.34 e B,-2- B fc>53»V > 



1/2 



^.exio 11 HzCe- 2 < 5 2 3 4(- 2 4 



(13) 
,(14) 



respectively. 

For v a < vx < v m , F vx oc t~i/ 4 . A steeper decline is 
possible (the steepest one is F Vx oc t~ 4 / 7 ), depending on the 
radiative correction, as shown in the upper panel of Fig. 2 
of Sari et al. (1998). 

The transition of the slow decline to a normal decline 
(F vx oc t~ 1/2 ) usually takes place at t ~ 0.1 day or earlier, 
when vx =v m . So we have 



C e ~ 0.016^^ 



^_ 2 t-^C P [2/(l + z)]^. 



(15) 



The numerical light curves is presented in Fig. 3. One 
can see that a long time multi-wavelength flattening is evi- 
dent with a small £ e - 

Before and after the temporal decline transition, the en- 
ergy spectrum of the XRT observation should be F v oc v~ x / 2 
and F v oc i/~ p / 2 , respectively. In other words, after the break 
in the light curve, the X-ray spectrum should be much softer 



3.3 Evolving shock parameters 

In the standard afterglow model, the shock parameters e e 
and es are assumed to be constant. However, it is also pos- 
sible that e e or e_g, or both, vary with time (see Yost et al. 
(2003) for detailed discussion). Fan et al. (2002) and Wei 
et al. (2006) modeled the optical flares detected in GRB 
990123 and GRB 050904 and found that both e e and e B 
of the forward shock (ultra-relativistic) and reverse shock 
(mild-relativistic to relativistic) were very different. This 
provides an indication evidence for a dependence of the 
shock parameters on the strength of the shock. Possible evi- 
dence for the shock strength dependent cb was also found by 
Zhang, Kobayashi & Meszaros (2003) , Kumar & Panaitescu 
(2003), McMahon, Kumar & Panaitescu (2004), Panaitescu 
& Kumar (2004), Fan, Zhang & Wei (2005b) and Blustin 
et al. (2006). Yost et al. (2003) and Ioka et al (2005) con- 
sidered afterglow emission assuming es and e e are time- 
dependent, respectively. Here we simply take (e e , es) oc 
(r~ a , r~ ft ) for T > r o , otherwise (e e , es) ~ const., where 
r o is the Lorentz factor of the outflow at the X-ray decline 
translation, both a and b are taken to be positive. For sim- 
plicity, we discuss only the case of a = b for T > r o . Below 
r o , the solution is the usual one. 

The typical synchrotron radiation frequency v m satisfies 

v m oc (r/r )- 5a/ V 3/2 oc t(i5- 24 >/« (i 6) 

where Y « 25^% [2/(1 + 3 \ 
The cooling frequency v c satisfies 



oc (r/r ) 3a/ V 1/2 oc t 



-(8+9a)/16 



(17) 
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Figure 4. The X-ray (0.2-10 kcV) and R-band afterglow light 
curves for the evolving shock parameter model. The parameters 
are listed in the figure. 



The maximum spectral flux F^, max satisfies 
F,, max oc(T/r o )- a/2 oct 3a/16 . 



(18) 



The observed flux behaves as (in this subsection, we take 
a = 1): 



i/jiiiax^c 



1/3 .*+?<• 
v c ' act 



for Vx < V c < V„ 



, 1/2 

iV.maxIV °C t 32 

-1/3 



for V c < VX < Vr, 



- e.max V m 



OC t 8 



0.38 



for v x <v m < v c \ 

P~ 1 15ap-9a-24p + 24 n „ 

F v ,m^v m 2 act ~t-°- 2 , 

for v m < v x < v c ; 

E=± ■ 16-18.-24p+15ap ^ t _ . 65 



1/2 
{ V n 



for fx > maxfi/c,!'™}. 



The afterglow light curves are shown in Fig. 4. As both 
e e and es increase with time, the flux of the early X-ray 
emission is dimmer than that of the constant shock param- 
eters model and the decline is much slower. Both are con- 
sistent with the current Swift XRT observations (Nousek et 
al. 2005). 



3.4 A very low nonconstant density 

In the standard ISM afterglow model, the number density of 
the medium is taken as a constant. In the wind model, the 
number density n decreases with the radius R as n oc R~ 2 
(Meszaros, Rees & Wijers 1998; Dai & Lu 1998; Chevalier 
& Li 2000). Here we discuss the general case n oc RT k (0 < 
k < 3). 

First, wc show that for a fireball decelerating in the 
BM self-similar regime (Blandford & McKee 1976), no X- 
ray flattening is expected regardless of the choice of k. The 
energy of the fireball is nearly constant and it is given by 
E iso w T 2 Mc 2 , where M oc R 3 ~ k is the total mass of 
the swept medium. So V oc i2~< 3 ~ fc )/ 2 . Considering that 



Figure 5. X-ray (0.2-10 keV) afterglow light curve for a very 
low nonconstant density: n 



10" 



(R/10 1 



" 3 for 10 16 < R < 10 19 cm and 



n = 10~ 7 cm -3 for R > 10 19 cm. 

dR oc F 2 dt oc R- (:i - k) dt, we have R oc tV(4-k) and 

r oc £-( 3 - fc )/[ 2 ( 4 - fc )]. 

Now v m decreases with t as 

v m oc r 4 ir fc / 2 oc r 3 ' 2 , (20) 

and v c and F^ jmax satisfy 

V C OC T- 4 R :ik/2 t- 2 OC t (3*-4)/P(4-*)] ; (21) 

F,, max oci? 3 - 3fe / 2 r 2 oct- fc /[ 2 < 4 - fc 'l, (22) 
j-rj9|pectively. This results in: 

t -30,-l)/4-fc/[2(4-fc)] for ^ m < ^ < ^ ; 



F 1 !/^ oc 



£-1/4 

t -(3p-2)/4 



for i/ c < < (23) 
for z^x > max{V m , ^ c }. 



The last two are independent of k. So no X-ray flattening 
appears. 

However, if the number density is sufficiently low, the 
deceleration timescale (oc n -1 ^ 3 ) can be very long and even 
as long as ~ 10 4 s. In this case, a slowly decaying X-ray 
afterglow may be obtained. One example has been plotted 
in Fig. 5, in which the density profile of the medium is taken 
as n = 10~ 4 cm" 3 for R < 10 16 cm, n = lO -4 ^ 1 cm" 3 for 
1 < Rie < 10 3 , and n = 10~ 7 cm" 3 for R 16 > 10 3 . An X-ray 
flattening appears when the shock front reaches R = 10 19 
cm. However, while the shape of the light curve is correct the 
X-ray flux is too low to account for most XRT light curves. 



3.5 Magnetized outflow 

A Poynting flux dominated outflow (Usov 1994; Thompson 
1994; Lyutikov & Blandford 2003) is an alternative to the 
standard baryonic fireball model. Within the context of this 
discussion it is of interest since it may also give rise to a 
slowly decaying X-ray afterglow (Zhang et al. 2005). We 
investigate, here, briefly this possibility, extended discussion 
will be presented elsewhere. 

We assume that the electromagnetic energy E p will be 



7 



transformed continuously into the kinetic energy of the for- 
ward shock. The dynamical evolution of the shocked medium 
is governed by (Huang et al. 2000; Wei et al. 2006): 



dF = - 



(F 2 - \)dm + dE p /c 2 



M oj + em + 2(1 - e)Fm ' ^ 

where E p = F 2 VB' 2 /(4n), V is the volume of the mag- 
netized outflow (measured by the observer) and B' is the 
comoving strength of the magnetic field. 

If the magnetic pressure is higher than the thermal 
pressure of the shocked medium, the magnetic pressure 
works on the shocked medium and the kinetic energy of 
the forward shock increases. A pressure balance between 
the shocked medium and the magnetized outflow is estab- 
lished, so we have (see also Lyutikov & Blandford 2003) 
B' 2 /(8n) — p gas ~ 4r 2 nm p c 2 /3, where p gas is the thermal 
pressure of the shocked medium. Therefore E p can be esti- 
mated by 4 

E p = 2F 2 P gas V « Sr 4 nmpC 2 V/3. (25) 

dEp/c 2 can be calculated as follows. Assuming that the 
whole system (the shocked medium and the magnetized out- 
flow) is adiabatic (i.e., the radiation efficiency e = 0), the 
energy conservation yields 

2F 2 p gas V + 3F 2 p gas (V to t -V) = Etot - r(M ej + m)c 2 , (26) 

where Vtot ~ 47ri? 2 A is the total volume of the system, A 
is the width of the system, which is described by dA = 
(P fsh -P)dR and /3 fsh ~ Vr 2 - 1/[F- is the velocity of 
the forward shock. Differentiating equation (26) we obtain: 

dEp/c 2 = 2{(16r 3 nmpVtot + M oj + m)dF + W^Rnrrip 



[R(P f3h 2A]dR + Fdm}. 



(27) 



After simple algebra, equation (24) can be rearranged 
as: (note that now we take e = 0) 



dF 



(r 2 + 2AT - l)dm + BdR 



(28) 



M oj + 2Tm + 2A(16T 3 nm p V t ot + M oj + m) ' 

where A = 1 for dE p /dR < 0, otherwise A = 0; B = 
32ATmm p r 4 R[(f3 fsh - (3)R + 2 A]. 

With proper boundary conditions and the relations 
dm = 4:TYnm p R 2 dR, dR = f3cF 2 (l + P)dt/(l + z), Vtot = 
47ri? 2 A and dA — ((3f s h — fi)dR, equation (28) can be solved 
numerically. In our numerical example, we take Ek = 10 52 
ergs, n = 1 cm~ 3 , E p — WEk and the width of the outflow 
is taken as 3 x 10 11 cm. The starting point of our calculation 
is at R — 2 x 10 16 cm (~ Rdec, the deceleration radius of 
the outflow, where ~ Ek/2 has been given to the shocked 
medium), at which F = 360 5 . We find out that most of the 
magnetic energy has been converted into the kinetic energy 

4 Providing that V oc R 2+t T- d , E p oc r 8 + 2c - d t 2 + c oc t' 5 (c, d, 
and 8 are all larger than 0), we have T oc £-(2+c+<S)/(8+2c-d) _ 
which should be flatter than r — (the canonical dynamical 
evolution of a ejecta without energy injection). It requires that 
2c + 3d < 8 — 8<5, otherwise dE p /dt < has been violated. It is 
evident that in the spreading phase, i.e., c = 1 and d = 2, E v 
can not be converted into the kinetic energy of the forward shock 
effectively. 

5 At that radius, the reverse shock has crossed the ejecta and a 
pressure balance between the shocked medium and the magne- 
tized outflow is reached. In this work, we do not calculate the re- 



of the forward shock in a very short time ~ 50(1 + z) s. 
A similar result has been obtained by Lyutikov & Bland- 
ford (2003). Though this timescale is much longer than the 
crossing time of the reverse shock, it is not long enough to 
account for the X-ray flattening detected in most GRBs. 



4 CASE STUDIES: CONSTRAINING THE 
MODELS 

GRB 050319 and GRB 050401, have well recorded X-ray 
and optical afterglows, with which the models discussed in 
§3 can be constrained. We discuss these constraints in de- 
tail here. For most Swift GRBs only the X-ray afterglow is 
well detected. Such bursts provide, of course, much weaker 
constraints on the model . We discuss one example, GRB 
050315, briefly. 



4.1 GRB 050319 

Both the optical and X-ray afterglows of GRB 050319 have 
been well recorded (Wozniak et al. 2005; Mason et al. 2005; 
Cusumano et al. 2005; Nousek et al. 2005). The optical flux 
declines with a power law slope of a = —0.57 between ~ 
200s after the burst onset until it fades below the sensitivity 
threshold of the UVOT after 5 x 10 4 s. The optical V-band 
emission lies on the extension of the X-ray spectrum, with an 
spectral slope /3 = —0.8 (Mason et al. 2005). The temporal 
behavior of the X-ray afterglow is more complicated. After 
a steep decay (a = —5.53) up to t = 370s, the light curve 
shows a slow decay with a temporal index of a = —0.54. 
It steepens to a — —1.14 at t = 2.60 x 10 4 s. The spectral 
indices in the slow decline phase and the normal decay phase 
are f3 = —0.7 and (3 = —0.8, respectively (Cusumano et al. 
2005; Nousek et al. 2005; However, see Quimby et al. 2006). 
Below we examine whether the models discussed above (in 
§3) can explain both the optical and the X-ray afterglows 
self-consistently. 

Energy injection: The energy injection model is be- 
lieved to able to explain the observation (e.g., Zhang et 
al. 2005; Mason et al. 2005; Cusumano et al. 2005). As 
shown in §3.1, for q — 0.6 and p = 2.4, both the optical 
and the X-ray afterglows decline as F vx oc t~ ' 54 when 
v m < vr < vx < Vci the corresponding spectral index 
should be (3 — — (p — l)/2 ~ —0.7. All these values are 
consistent with the observation. However, the non-detection 
of the further X-ray break caused by the spectral transla- 
tion (y c < vx) up to ~ 10 s after the trigger suggests that 
e B ~5x 10~ 3 and n ~ 10~ 3 cm~ 3 (Cusumano et al. 2005). 
The problem is that F Vx depends on n and e_g sensitively 
for vx < v c (see equation (7)). The smaller n and es, the 



verse shock emission (see Fan, Wei & Wang (2004a) for the reverse 
shock emission with mild magnetization and Zhang & Kobayashi 
(2005) for reverse shock emission with arbitrary magnetization). 
With the ideal MHD jump condition, the reverse shock can not 
convert the magnetic energy into the kinetic energy of the for- 
ward shock effectively, as shown in Kennel & Coronitti (1984), 
Fan, Wei & Zhang (2004b) and Zhang & Kobayashi (2005) both 
analytically and numerically. 
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smaller F vx . We show below that it is quite difficult to re- 
produce the detected X-ray and optical light curves with the 
energy injection model. 

The earliest R-band data is collected at ~ 200 s (note 
that the real onset of GRB 050319 is about 130 s before 
the Swift trigger time taken in Woznik et al. (2005), see 
Cusumano et al. (2005) for clarification), and the flux is 
about F VR ~ 0.7 mjy. At that time, the total energy of the 
outflow is still dominated by the initial Ek , and -F„, max and 
v m are still described by equations (2) and (3) , respectively. 
The conditions F„ ilnax > 0.7 mjy and v m (t ~ 200 s) < vr 
yield 

e^_ 2 E k , bZ n\l 2 > 0.8 => E k , 53 > 0.8e-^ 2 2 n" 1/2 , (29) 
<5 3 4 /2 - 2 4-i < 0-04 => e e <0.02E^i 4 e B % (30) 

respectively. The condition v c > vx ~ 10 17 Hz holding up 
to t ~ 10 6 s gives 

EkgeB%W{l + Y )- a >94Q 

^e B < W- 4 E^n- 2/ *(l + Y o y 4 / 3 , (31) 

where Y is the Compton parameter at t ~ 10 6 s. To derived 
this relation, we assume that at t ~ 200 s, v c is described by 
equation (4), and v c oc t (q ~ 2)/2 ~ i~ ' 7 up to t ~ 2.6 x 10 4 
s (i.e., in the energy injection phase), then v c oc i -1 / 2 up to 
t ~ 10 6 s. Combing equations (29-31), we have 

£fc,53 > 13«o 1/5 (l + Y ) 4/5 , (32) 
e e < 0.06£-^ 6 nJ /6 (l + y o ) 1/3 . (33) 

Now Yo ~ (see the Appendix for detail), we have 
Sfe > 1.3 x 10 54 ergs n^ 1 ^ . On the other hand, the en- 
ergy injection coefficient A' = Ac 2 ~ (1 + z)Ek/to ~ 
1.4 x 10 S2 rio 1/8 ergs s" 1 for t ~ 370 s. Please note 
that .A' is comparable to the recorded luminosity of most 
GRBs and the X-ray luminosity recorded by XRT is just 
~ 10 48 ergs s _1 . The outflow accounting for the late time 
injection is so energetic that strong soft X-ray to 7— ray 
emission powering by shocks or magnetic dissipation are ex- 
pected. They will quite likely dominate over the correspond- 
ing forward shock emission, which is inconsistent with the 
observation. 

This model is also disfavored by the different temporal 
behavior of the X-ray and the optical afterglows at t > 2.6 x 
10 4 s. We therefore conclude that the energy injection model 
can't account for the multi-wavelength afterglows of GRB 
050319. We tried to fit both the R-band and X-ray afterglows 
with reasonable parameters numerically but failed. 

Provided that the energy injection model works (i.e., 
there is a mechanism to keep such energetic outflow steady 
enough and there is no magnetic dissipation), the initial 
GRB efficiency in this case is as low as e 7 = E 1 / [E 1 + Ek) ~ 
0.08 nl /5 . 

Small f e : This model is disfavored by two facts. One is 
that in the X-ray flattening phase, v c < v x < v m , the corre- 
sponding spectral index is /3 ~ —1/2, which only marginally 
matches the observation ~ —0.7. The other is that after the 
temporal transition at t ~ 10 4 s, i/ c < % < vx, the spec- 
tral index should be /3 — — p/2 ~ —1.2, which is inconsistent 
with the observation. 

Evolving shock parameters: As shown in §3.3, for 
v m < vr <vx < v c , p = 2.4 and a = b — 0.6, (F VR , F vx ) oc 



t~ ' 54 and the spectral index [3 = -{p- l)/2 0.7, are all 

consistent with the observation. After the shock parameters 
saturate at t ~ 2.6 x 10 4 s, F Vx oc t~ 1A and f3 = —0.7 as 
long as vx < v c , which also matches the observation. How- 
ever, the optical light curve should be much steeper since 
v m < vr < Vc also holds. The UVOT observation and the 
ground based R-band observation suggest that the decline of 
optical emission does not change up to t ~ 2 x 10 5 s, though 
the scatter of the flux is quite large (Kiziloglu et al. 2005; 
Sharapov et al. 2005; see Mason et al. 2005 for a summary). 
Therefore the evolving shock parameter model is disfavored. 

Very low nonconstant density: With proper param- 
eters as well as proper density profile, an X-ray flattening 
does appear (see Fig. 5). However, as already mentioned, the 
flux is too low to match most observations, here we do not 
discuss it further. 

Off-beam annular jet model. Recently, Eichler & 
Granot (2005) suggested that the flat part of the XRT light 
curve may be a combination of the decaying tail of the 
prompt 7— ray emission and the delayed onset of the after- 
glow emission observed from viewing angles slightly outside 
of the edge of the jet (i.e., off-beam). This model, like oth- 
ers mentioned above, can account for the slow decline of 
many X-ray afterglows, but may be unable to explain both 
the optical and the X-ray afterglows of GRB 050319 self- 
consistently, as shown below. 

Following Eichler & Granot (2005) , we assume that the 
off-beam angle is 86 ~ l/rj n t, where Ti nt is the initial 
Lorentz factor of the outflow. Larger 86 is less favored since 
the slowly decayingR-band afterglow has been well recorded 
as early as t ~ 200 s, which implies that the afterglow on- 
set has not been delayed so much. In the off-beam case, the 
typical synchrotron radiation frequency should be 

,„„ 1n ll u u.1/2 1/2 2 ^-,2/ 1 + z O/2 1/2,-3/2 , OA s 

v m ^ 7.6x10 HzS k ; 53 e B / _ 2 e ei _iC p (— — ) ' a 1 t d ,(34) 

where a [1 + {V int 86) 2 ] ~ 2 is the Doppler factor. There- 
fore the condition v m (t ~ 200 s) < vr results in 

U < O.OnE-^e-^ia^y 1 ^. (35) 

For 86 ~ 1/Ti n t, the late time (i.e., the normal decline 
phase) afterglow emission is quite similar to the on-beam 
case (Eichler & Granot 2005). 

We use equation (7) to estimate the late time X-ray flux, 
though the predicted flux of the annular jet model should 
be somewhat different from that of our conical jet model 
(Granot 2005; Eichler & Granot 2005). The XRT flux ^ 
8 x 10~ 12 ergs s _1 cm~ 2 at td ~ 0.3 gives 

£ fe , 53 « 0.33e 4( i7 rt/(j,+3) e s ^+ 1)/(p+3) n - 2/(j,+3) . (36) 

The condition v c > v x ~ 10 17 Hz holding up to t ~ 10 6 s 
yields 

e s <3xlO- 4 £ fe l/ 3 n- 2/3 . (37) 

Combing equations (35-37), we 

have E k > 10 55 ergs n ( ; 1/5 (a/2) 3(p " 1)/10 . While we man- 
age to fit both X-ray and optical data, the energy needed is 
too large for any realistic progenitor models. We therefore 
suggest that the off-beam annular jet model is also unable 
to account for the afterglows of GRB 050319. 
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4.2 GRB 050401 

The early X-ray light curve is consistent with a broken power 
law with a = —0.63 and —1.41 respectively, the break is 
at t b ~ 4480 s (de Pasquale et al. 2006). The X-ray spec- 
tral indices before and after the break are nearly constant 
~ —0.90. Therefore the small £ e model is ruled out directly. 
Zhang et al. (2005) also show that the flat electron distribu- 
tion model (1 < p < 2) is unable to account for the X-ray 
afterglow observation. The afterglow has also been detected 
in R-band, which decays as a simple power law oc t~ ' 76 up 
to i~3.5x 10 4 s (Rykoff et al. 2005). 

Energy injection (p ~ 2.8). After the break, the light 
curve is consistent with an ISM model for v m < fx < v c 
with p ~ 2.8. Before the break, it is consistent with the same 
model with q — 0.5 (see also Zhang et al. 2005). As far as 
the R-band afterglow emission is concerned, there are two 
possibilities. One is that v m < vr < v c , the optical afterglow 
should follow the temporal behavior of the X-ray afterglow, 
which is not the case. The other is that vr < v m for t < tb, 
the afterglow increases as t for q ~ 0.5 (see §3.1), which 
is inconsistent with the observation. We therefore conclude 
that the popular energy injection model is unable to account 
for the data in this burst as well. 

Evolving shock parameters (p ~ 2.8). The light 
curve after the break is consistent with an ISM model for 
Vm <vx <v c with p ~ 2.8. Before the break, it is consistent 
with the same model with a = b — 0.7. Can it reproduce 
the optical afterglow? The answer is negative. Provided that 
vr < v m for t < tb, the optical afterglow should increase as 
t 0A (see §3.3), which is inconsistent with the data. The case 
of v m < vr is ruled out directly in view of the different 
temporal behavior of X-ray and R-band afterglows. 



4.3 GRB 050315 

After a steep decay up to tbi = 308 s, the X-ray light curve 
shows a flat "plateau" with a temporal index of a = —0.06 
(the spectral index of XRT data is f3 = —0.73). It then 
turns to a — —0.71 at tb2 = 1.2 x 10 4 s, the spectral index is 
13 = -0.79. Finally there is a third break at t b3 = 2.5 x 10 5 
s, after which the temporal decay index is a = —2.0 and the 
spectral index is (3 = —0.7 (Nousek et al. 2005; Barthelmy 
et al. 2005). 

There are two possible interpretations for the long 
term constant spectral index (3 ~ —0.7. One is that 
mayi{v c ,v m } < vx after t = 308 s and the power law in- 
dex of the shocked electron p ~ 1.5. The other is that 
v m < vx < v c for tbi < t < tb3 and p ~ 2.5. 

Energy injection (p ~ 2.5). To obtain the slow decline 
for t b i < t < hi, energy injection with q ~ 0.2 is needed, q ~ 
0.9 is needed to reproduce the X-ray afterglows at t b i < t < 
tbz- The late time sharp decay appears when the boundary 
of a non-spreading jet becomes visible. 

Evolving shock parameters (p ~ 2.5). As shown in 
§3.3, with a — b — 1.2, we have a slow decline slope a = 
—0.06 between tbi and tbi- To get a decline slope a = —0.71 
between tbi and t b z, a — b = 0.45 are needed. The late time 
sharp decay appears when the boundary of a non-spreading 
jet becomes visible and a = 6 = 0.45. 

We find that both models can explain the observed X- 
ray light curves of GRB 050315. 



5 SUMMARY & DISCUSSION 

During the past several months, the Swift XRT has collected 
a rich sample of early X-ray afterglow data. A good frac- 
tion of these afterglows show a slow decline phase lasting 
between a few hundred to several thousand seconds. The 
energy injection model is the leading model to account for 
these slowly decaying X-ray afterglows (e.g., Zhang et al. 
2005; Nousek et al. 2005; Panaitescu et al. 2006; Granot & 
Kumar 2006). It has been suggested that in this model, the 
GRB efficiency might be as high as 90%. Such a high effi- 
ciency challenges the standard internal shock model for the 
prompt 7— ray emission. 

In this work, we have re-examined the GRB efficiency of 
several pre-Swift GRBs and one Swift GRB. In addition, we 
have explored several mechanism which might give rise to 
a slowly decaying X-ray light curve and we have compared 
the predictions of these models with the well recorded multi- 
wavelength afterglows of GRB 050319 and GRB 050401. We 
draw the following conclusions: 

1. The GRB efficiency of pie-Swift GRBs that has been 
derived directly from the X-ray flux 10 hours after the burst 
has been overestimated. For these Swift GRBs with long 
time X-ray flattening, the GRB efficiency is also moderate 
(around 0.5), even when taking into account the possibility 
of energy injection. Such efficiency can be understood within 
the standard internal shock model. 

2. With a proper choice of parameters, the slow decline 
slope of X-ray afterglow like the one detected in GRB 050319 
can be well reproduced by several models — the energy injec- 
tion model, evolving shock parameter model (in which the 
shock parameters are assumed to increase with the decrease 
of the shock strength for t < 10 4 s), the small £ e model 
(in which the shock energy has been give to a fraction £ e 
of electrons, rather than total) and the very low noncon- 
stant density model. Out of these models, the last two are 
ruled out by the X-ray data itself. In the last model, the 
resulting X-ray afterglow is too dim to match most XRT 
observations. The small £ e model is also disfavored since (1) 
In the slow decline phase, the XRT spectrum are usually 
much softer than v~ x ^\ (2) After the light curve break, no 
spectral steepening has been detected in most cases, which 
is inconsistent with the model. The other models, including 
the energy injection model and the evolving shock param- 
eter model seem to be consistent with the X-ray afterglow 
observations. 

3. While two models: the energy injection model and 
the evolving shock parameter model are consistent with the 
X-ray data, they fail to reproduce both the X-ray and the 
optical afterglows of GRB 050319 and GRB 050401. In each 
burst, the optical flux declines slowly up to ~ 10° s. On 
the other hand, the X-ray light curve decays slowly up to 
t ~ 10 4 and then turns to the normal faster decay (F oc i~ 1,2 
or so). The temporal index of the slow decay X-ray phase is 
close to that of the optical light curve. The XRT spectrum 
is unchanged before and after the X-ray break. This means 
that the break is not caused by a cooling break in which v c 
crosses the observed frequency. . 

The failure of all models that we considered to fit both 
the X-ray and the optical afterglow light curves suggests 
that we should look for another alternative. An intriguing 
possibility is based on fact that the extrapolation backwards 
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of the late X-ray light curve is in agreement with (or 1 to 2 
order lower than) the prompt X-ray emission. This suggests 
that we face a "missing energy problem". Namely, during 
the slow decay phase (in which the X-ray flux is rather low) 
we miss X-ray emission. Is it possible that during this phase 
this energy is dissipated into a different channel and not 
into Synchrotron X-rays and that this different channel be- 
comes ineffective at around ten hours? Put differently, dur- 
ing this phase the electrons within the forward shock emit 
Synchrotron X-rays inefficiently. A possibility of this kind 
(that we have considered and found not to work) is if the 
X-ray emitting electrons are cooled efficiently via inverse 
Compton (and hence their Synchrotron X-ray emission is 
weaker). As already mentioned inverse Compton cooling is 
important in determining the X-ray flux. Furthermore, due 
to the Klein-Nishina cutoff this cooling becomes unimpor- 
tant at approximately one day. However, this transition is 
not sharp enough to produce the observed slowly decaying 
X-ray light curves. While inverse Compton cooling does not 
work it is possible that another, yet unexplored, process of 
this kind is responsible for the observed light curves. 



ACKNOWLEDGMENTS 

We thank A. Panaitescu and E. Waxman for their construc- 
tive comments. We also thank the referee for helpful sugges- 
tion. Y. Z. Fan thanks D. M. Wei for his help and J. Dyks for 
checking Fig. 5. This work is supported by US-Israel BSF. 
TP acknowledges the support of Schwartzmann University 
Chair. YZF is also supported by the National Natural Sci- 
ence Foundation (grants 10225314 and 10233010) of China, 
and the National 973 Project on Fundamental Researches of 
China (NKBRSF G19990754). 



REFERENCES 

Akerlof C, ct al. 1999, Nature, 398, 400 
Barthelmy S. D., ct al. 2005, ApJ,635, L133 
Blandford R. D., & McKee C. F. 1976, Phys. Fluids., 19, 1130 
Beloborodov A. M., 2000, ApJ, 539, L25 
Berger E., Kulkarni S. R., Frail D. A., 2004, ApJ, 612, 966 
Blustin A. J., et al. 2006, ApJ, 637, 901 
Campana S. et al. 2005, ApJ, 625, L23 
Chevalier R. A., Li Z. Y., 2000, ApJ, 536, 195 
Cusumano G., ct al. 2005, ApJ, in press (astro-ph/0509689) 
Dai Z. G., 2004, ApJ, 606, 1000 
Dai Z. G., Lu T., 1998a, A&A, 333, L87 
Dai Z. G., Lu T., 1998b, MNRAS, 298, 87 
Daigne F., Mochkovitch R., 1998, MNRAS, 296, 275 
de Pasquale M. et al. 2006, MNRAS, 365, 1031 
Dyks J., Zhang B., Fan Y. Z., 2005, ApJ, submitted (astro- 
ph/0511699) 

Eichler D., Granot J., 2005, ApJ, submitted (astro-ph/0509857) 
Fan Y. Z., Dai Z. G., Huang Y. F., Lu T., 2002, Chinese J. Astron. 

Astrophys., 2, 449 (astro-ph/0306024) 
Fan Y. Z., Wei D. M., Wang C. F., 2004a, A&A, 424, 477 
Fan Y. Z., Wei D. M., Zhang B., 2004b, MNRAS, 354, 1031 
Fan Y. Z., Zhang B., Proga D., 2005a, ApJ, 635, L129 
Fan Y. Z., Zhang B., Wei D. M., 2005b, ApJ, 628, L25 
Freedman D. L., Waxman E. 2001, ApJ, 547, 922 
Granot J. 2005, ApJ, 631, 1022 

Granot J., Konigl A., Piran T., 2006, MNRAS submitted (astro- 
ph/0601056) 



Granot J., Piran T., Sari R., 1999, ApJ, 513, 679 
Granot J., Kumar P., 2006, MNRAS, 366, L13 
Granot J., Nakar E., Piran T., 2003, Nature, 426, 138 
Huang Y. F., Gou L. J., Dai Z. G., Lu T., 2000, ApJ, 543, 90 
Ioka K., Toma K., Yamazaki R., Nakamura T. 2005, A&A sub- 
mitted (astro-ph/0511749) 
Kennel C. F., & Coronitti E. V., 1984, ApJ, 283, 694 
Kiziloglu U., et al., 2005, GCN Circ. 3139 
Kobayashi S., Piran T., Sari R., 1997, ApJ, 490, 92 
Kobayashi S., Sari R., 2001, ApJ, 551, 943 
Kumar P. 2000, ApJ, 538, L125 
Kumar P., Panaitescu A., 2003, MNRAS, 346, 905 
Kumar P., Piran T., 2000, ApJ, 532, 286 
Llod-Ronning N. M., Zhang B., 2004, ApJ, 613, 477 (LZ04) 
Lyutikov M., Blandford R., 2003 (astro-ph/0312347) 
McMahon E., Kumar P., Panaitescu A., 2004, MNRAS, 354, 915 
Mason K. O., ct al., 2005, ApJ, in press (astro-ph/0511132) 
Meszaros P., Rees M. J., Wijers R. A. M. J., 1998, ApJ, 499, 301 
Modcrski R., Sikora m., Bulik T., 2000, ApJ, 529, 151 
Nousek J. A., et al., 2005, ApJ, in press (astro-ph/0508332) 
Paczynski B., Xu G. H., 1994, ApJ, 427, 708 
Panaitescu A., Kumar P., 2002, ApJ, 571, 779 
Panaitescu A., Kumar P., 2004, MNRAS, 353, 511 
Panaitescu A., Meszaros P., Gehrels N., Burrows D., Nousek J. 

2006, MNRAS, in press (astro-ph/0508340) 
Panaitescu A., Meszaros P., Rees M. J., 1998, ApJ, 503, 314 
Papathanassiou H., Meszaros P., 1996, ApJ, 471, L91 
Piran T., 1999, Phys. Rep., 314, 575 
Quimby R. M., et al., 2006, ApJ, 638, 920 
Rykoff E. S., et al. 2005, ApJ, 631, L121 
Rees M. J., Meszaros P., 1994, ApJ, 430, L93 
Rees M. J., Meszaros P., 1998, ApJ, 496, LI 
Sari R. 1997, ApJ, 497, L17 
Sari R., Meszaros P., 2000, ApJ, 535, L33 
Sari R., Narayan R., Piran T., 1996, ApJ, 473, 204 
Sari R., Piran T., 1997a, MNRAS, 287, 110 
Sari R., Piran T., 1997b, ApJ, 485, 270 
Sari R., Piran T., Narayan R. 1998, ApJ, 497, L17 
Sharapov D., Ibrahimov M., Karimov R. , Kahharov B., Poza- 
nenko A., Rumyantsev V., Beskin G., 2005, GCN Circ. 3140 
Thompson C, 1994, MNRAS, 270, 480 
Usov V. V., 1994, MNRAS, 267, 1035 
Vaughan S. et al. 2005, ApJ, in press (astro-ph/0510677) 
Wei D. M., Lu T., 1998, ApJ, 505, 252 
Wei D. M., Lu T., 2000, A&A, 360, L13 
Wei D. M., Yan T., Fan Y. Z., 2006, ApJ, 636, L29 
Wijers R. A. M. J., & Galama T. J. 1999, ApJ, 523, 177 
Wozniak P. R., Vestrand W. T., Wern J. A., White R. R., Evans, 

S. M., Casperson D., 2005, ApJ, 627, L13 
Yost S., Harrison F. A., Sari R., Frail D. A., 2003, ApJ, 597, 459 
Zhang B., Fan Y. Z., Dyks J., Kobayashi S., Meszaros P., Burrows 
D. N., Nousek J. A., Gehrels N. 2005, ApJ, in press (astro- 
p/0508321). 
Zhang B., Kobayashi S., 2005, ApJ, 628, 315 
Zhang B., Kobayashi S., Meszaros P., 2003, ApJ, 595, 950 
Zhang B., Meszaros P., 2001, ApJ, 552, L35 
Zhang B., Meszaros P., 2002, ApJ, 566, 712 



APPENDIX A: THE GENERAL FORM OF THE 
INVERSE COMPTON PARAMETER 

For the photons with frequency higher than v, the Compton 
parameter should be suppressed significantly since it is the 
Klein-Nishina regime, where v is governed by (1 + z)~jj\v ~ 
rm e c 2 , i.e., 

z>~ 1.2 x 10 20 Hz (l + z) -1 !^ 1 . (Al) 
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We extend the derivation of the Compton parameter Y 
given by Sari et al. (1996) to the general form, in the limit 
of single scattering. The ratio of the inverse Compton power 
(Pic) to the synchrotron power (P ay n) of an electron with 
random Lorentz factor 7 e is given by 

K(7e) ~P~- ~U^~ - [1 + F( 7e )]e B ' [A2) 

where rj KN is the fraction of synchrotron radiation energy of 
total electrons emitted at frequencies below v. So we have 



F( 7e ) = (-1 + V 1 + 4 WjfJv e e /€ S )/2. (A3) 
Below we estimate the parameter r\ KN in different cooling 



regimes. 



A. Slow cooling. 

(z,K)- (p - 1)/2 , 



for v m < v < v c ; 
for v c < v < vm- 



(A4) 



where vm ~ 2.8 x 10 22 r/(l + z) Hz is the maximal syn- 
chrotron radiation frequency of the electrons accelerated by 
the forward shock. For p > 2, the total energy emitted 



(p-l)/2r 



,(3-p)/2 



,(3-p)/2 1 



where 



the photons with frequencies below v m have been ignored. 
Throughout the Appendix, v c is still described by equation 
(4) but without the correction of 1/(1 + Y) 2 . We have 



0, 



g (3 -p)/2_ 1/ (3-p)/2 

„(3-p)/2™ (3-p)/2 1 



, _ (3-p)^ 1/2 g( 2 -P)/ 2 
1 (3-p) (3-P) 1 

^ c 2 -(p-2)„ m 2 



for £ < v m ; 

for ^ m < v < !/ c 

for < v. 



(A5) 



For 1 < p < 2, the total energy emitted is J F v dv 
Si, where Si = [(3 - p)^ c 1/2 ^ p)/2 



g_Fb 



>-l)/2 



(2-p)(3-p)" c 

M (3-p)/2 _ ^ _ p)i/£? _p ^ 2 ]. Now r/ KN can be estimated as 



f 0, 

(2-p)(f.< 3 -P)/ 2 -^- p)/2 ) 
Si 

1 (2-p) (2-p) 

, _ CS-v)v?{y M ' 1 -0 2 ) 

4 C. 1 



for £> < v m ; 

for v m <v < v c ; 

for i^c < z> < i^m- 



APPENDIX B: WHEN IS THE KLEIN-NISHINA 
CORRECTION IMPORTANT? 

In the shock front, the magnetic field strength B is 

B = 0.04r4 /2 _ 2 nJ /2 . (Bl) 

The typical synchrotron radiation frequency of an electron 
with random Lorentz factor y e is 

v( le ) = 2 -^-Bz^rB, (B2) 

1 + 2 

A. The XRT lightcurve 

For v x ~ 10 17 Hz, we have ^(yx) = 1.3 x 10 5 [2/(1 + 

^"'^rr^B^W 1/4 and 



?~5x 10 15 Hz [(l + z)/2] 1/2 r?e 



l/2 F 2 £ l/4 1/4 
B,-2"0 ' 



(B3) 



Therefore, the Klein-Nishina correction seems to be 
unimportant (i.e., ri KN ~ 1) for t<2 ~ 1 (when z> ~ f c ) and 
e s ~0.01. 

B. The R-band lightcurve 

For v R ~ 4.3 x 10 14 Hz, we have 7 e (^ii) = 8 x 10 3 [2/(1 + 



-l/2 r -l f 



-1/4 -1/4 



and 



8x 10 la Hz [(l + z)/2] 1/2 r^e 



l/2 p 2 £ l/4 1/4 
-2 ,t 



(B4) 



Then, with es ~ 0.01, the Klein-Nishina correction 
seems to be unimportant for a long time. On the other hand, 
the factor r] ~ min{l, (^ m /^) (p ~ 1)/2 } ~ 1 for t d < 1. As a 
consequence, the inverse Compton effect is very important 
both for the long wavelength afterglow calculation and for 
the X-ray lightcurve calculation. However, it may be unim- 
portant for a lower es since v c oc e s 3 ^ 2 . 



(A6) 



B. Fast cooling. 

. , , W^)- 1/2 , 



for v c < v < v m ; 
for v m < v < vm- 



(A7) 



For p > 2, the total energy emitted is J F v dv = 

2Foi/m[(jE^)vm 2 — v\ ], where the emission below v c has 
been ignored. The r/ KN is estimated as 

0, for v < v c ; 

for v c < v < v„ 



»!/J_ 1/2 



(f^)^ /2 



1/2 1 
(p-l)/2 ( 2 _p)/2 



(A8) 



7 i \ !/ 2 "{ n\ !/ 2 ' 



for ^ m < v. 



For 1 < p < 2, the total energy emitted is J F v dv 



1/2 



2-P 
1/2, 



-S2, where S2 = [v, 
We have 

r 0, 

(2-p)( i >V 2 _„V 2 ) 



(p-l)/2 (2-p)/2 



(p-l)^ /2 -(2- 
for v < v c ; 



